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Large Mass and Chemical Potential Model: A Laboratory for QCD? 

Ralf Hofmann^and lon-Olimpiu Stamatescu"'*^ 

^Institut fiir Theoretische Physik, Univ. Heidelberg, Heidelberg, Germany 
^Forschungsstatte der Evang. Studiengemeinschaft, Heidelberg 

We use a model based on the hopping parameter expansion to study QCD at large fi. We find interesting 
behavior in the region expected to show flavor-color locking. 



Motivation and problem setting While the 
analysis for QCD at small /i and large T gradu- 
ally improves, there are no studies yet at large fi 
and small T. Results hereto have only been ob- 
tained in models which are not directly derivable 
from QCD, such as SU(2) and isovector chemical 
potential, or NJL. Here we want to approach this 
very interesting region ''T] in a model which can be 
obtained as large mass, large chemical potential 
approximation of full QCD |2| and which in the 
static, dense quark limit T has been proposed 
as quenched formulation of QCD at finite baryon 
density See Fig. For this we expand the 
fermion determinant in the inverse mass to sec- 
ond order, which ensures a limited mobility of the 
quarks. We use Wilson fermionic action and have 
for the hopping parameter: 

K= [2{M + 3 + cosh fi)]-^ = [2(Afo + 4)]"i 

(M : bare mass, Mq : bare mass at /i = 0). Then: 



zP{K,ti,{U}) 



exp 




r,q,i.t,t' 

^ ZP {C, {[/}) n Det (l - (e Cy V^J,^, 

x,r,q,i,t,t' 



(1) 



(2) 



In the temporal gauge (Un,4 — 1, except for 
U{s,7i4=N^}A = ^s- f^'se) we have P3 = V3 and 

-P^Lt' = iVsr-'U^,,t)AVs+,rU^,,r), (3) 



with r > q >0, i = ±1,±2,±3, 1 < t < t' < Nr 
{t < t' for q = 0). Here e = 1 (-1) for p.b.c. 
(a.p.b.c). For detail see 

The hope is that the model retains some fea- 
tures of QCD in the large ^ region, and that due 
to the simplicity of the determinant large statis- 
tics reweighting methods could converge, giving 
insight in the phase structure of QCD. 



Fixed mass plaj 



\ 7 














/ 








1 













SniaJl 
QCD 



I A ; Quenched QCD ; 
1^ ~ 250 MeV 




_ _ 



7"=tt fuU QCD 



K - 1 / Ma; 



Figure 1. Tentative phase diagram of QCD. 



One flavor QCD The simulations used here and 
below are done on a G'* periodic lattice (a.p.b.c. 
for fermions in r-direction, Wilson plaquette ac- 
tion for the gauge field). The algorithm uses a 
local Boltzmann factor obtained by taking out 



2 



from Zp the factor 



most generally is parameterized as 



JJexp |2C7^eTr 



r.q.i,t,t' 



■pr,q 

X'L.t.t' 



compensating for this in the global reweighting 
(both procedures are vectorized). We use tem- 
poral gauge fixing for easy book-keeping of the 

[21 

contributions to Zp . We first analyze one fiavor 
QCD, some results are shown in Fig. [3 The con- 
vergence appears rather good, up to very large ns 
(baryon density over T^; nsisat) — 216). Both 
ub and Polyakov loop increase above 1. 
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Figure 2. One flavor QCD. Upper plot: itera- 
tion history for n^. Lower plot: ns /nB{sat) and 
Polyakov loop vs ^. 

Application to QCD with three degener- 
ate flavors: Color-flavor locking We now in- 
vestigate whether a transition to a phase, where 
the ground state is characterized by a color-flavor 
locking condensate of quark bilinears PP, takes 
place. Disregarding the possibility of parity vi- 
olation in the ground state, such a condensate 



K2 s^st 



(4) 



where C = ^7270 is the charge conjugation ma- 
trix. The contraction of Dirac indices is implicit 
while a ■ ■ ■ h and i ■ ■ ■ z denote color and flavor 
indices, respectively. The condensate is invari- 
ant under locked global color and vectorial fla- 
vor rotations. To project out the invariants Ki 
and K2 we apply \ [SfS^ij + pS^S^) to the con- 
densate (01) with (A,p) — (—^,—3) and (A,p) = 
(i,— i), respectively. We define the two-point 
correlator of color-flavor locking quark bilinears 



[^PU^Wi^ji^) [rkiO)Cl'^tm ')■ Two- 
point correlations of a certain combination of the 
invariants K12 are obtained by applying 



A(x)A(0) {6f6'^+p{x)5^6'^) {SlSf 



(5) 



to the above correlator after accordingly adjust- 
ing (A(a;), p{x)) and (A(0), p(0)). In the Euclidean 
formulation the result can be decomposed as 

A(x)A(0) J2 {i^n^)Oif''i^',ix)^t{0)Oifrk{0) 

I, .1=1,2 

= O'^BO-ho [w^^{x,Q)Wb^d{^,^) 

- WA^j,{x,Q)WB^ci^,Q)) deiW (6) 



(see, e.g., see for example 0), where 

= <5r 5,^27475, Olf = p{x)5';5h2lAl^ , 



V 

0~l,cd seed r\2,c 

ki = (5^(5,757472, 0^:^ 



Here a 'super' index A = (a, a, i) was introduced 
with a a Dirac index. Like det W W"^ can also 
be expanded in powers of k [2], the paths con- 
tributing to C(k^) are shown in Fig. |3| 

We have performed simulations for QCD with 
3 degenerate flavors at /? = 5.6 on a 6^* lattice 
at K = 0.12 and various p (same for all flavors). 
Results are shown in Fig. ^ The general con- 
vergence is poorer than for the 1-flavor case at 
same p. In particular, the region around /i = 1 
acknowledges strong fluctuations, while at larger 
p the situation improves again. Here we also 
measure the "QQ" -susceptibility obtained by in- 
tegrating the correlator of K\ + ^^^2 ©, where 
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Contributions to the quark propagator to order K C 
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Contributions to the di-quark propagator to order K ^ 

Figure 3. Paths contributing to 0{k'^). 



for definiteness we chose ^ = 0.5. For this quan- 
tity we use maximal gauge fixing. The contri- 
bution appears essential for the non-trivial large 
fi behavior (the "mobility" defined as the 0{k^) 
fraction of ub over total ns is typically ^ 0.3). 

Conclusions and outlook We have performed 
simulations at small T, large fi for a model which 
can be obtained from QCD by small order hop- 
ping parameter expansion at large /i. The quarks, 
although very heavy, have some limited amount of 
mobility. The results show strong increase of the 
baryon density and other observables above ^ 1 
which may be a signal for entering a new, high 
density phase at small temperature. Particularly 
interesting is the behavior of the di-quark correla- 
tor in 3-flavor QCD which becomes increasingly 
flat at large /z, leading to a strongly increasing 
susceptibility. This may be a signal for the devel- 
opment of a condensate with color-flavor locking. 
Simulations on larger lattices would be essential 
to check this conjecture, it remains to be seen 
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Figure 4. 3-flavor QCD. Upper plot: Iteration 
history for at = 0.85 (multiplied by 10) and 
1.08. Lower plot: "QQ" -susceptibility to 0{k°) 
(multiplied by 5) and to O(k^), vs fj,. 

whether good convergence can be achieved in that 
case. 
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